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We present a new technique to extract information on the Unitarity Triangle from the study of 
B — > Knn Dalitz plots. Using the sensitivity of Dalitz analyses to the absolute values and the 
phases of decay amplitudes and isospin symmetry, we obtain a new constraint on the elements of 
the CKM matrix. We discuss in detail the role of electroweak penguins and outline future prospects. 



INTRODUCTION 

The study of flavour and CP violation in B decays 
allows us to test the flavour structure of the Standard 
Model (SM) and to look for New Physics (NP). In the 
last few years, B factories have provided us with a large 
amount of new data in this field, in particular the first 
measurements of the angles a and 7 of the Unitarity 
Triangle (UT). These measurements are in agreement 
with the indirect determinations of the UT angles Q, Q 
and therefore provide an important test of the Standard 
Model and stringent constraints on NP 0. 

B — > Kir decays could in principle constitute an- 
other source of information on 7 jj, |5|. However, the 
fact that the tree-level b — > u transition, carrying the 
phase 7, is doubly Cabibbo suppressed in these chan- 
nels, together with a large dynamical enhancement of the 
Cabibbo allowed penguin contribution Q , make a model- 
independent extraction of 7 from B — » Kn impossible. 
Model-dependent studies of these channels have tackled 
this issue with no success 0, • 

In this letter, we consider the possibility of obtaining 
a tree-level determination of 7 from B — > K*ir decays, 
cancelling the effect of penguins through the rich set of 
information available from Dalitz plot analyses. In par- 
ticular, the B factories can provide the magnitude and 
phase of decay amplitudes separately for B and B de- 
cays. One can think of exploiting this information, to- 
gether with isospin symmetry, to build combinations of 
amplitudes that are proportional to a single weak phase. 

This very simple idea, however, in the case of B — > 
Ktttt, has to face the presence of Electroweak Pen- 
guins (EWP), doubly Cabibbo enhanced with respect to 
current-current operators in b — > s transitions. This en- 
hancement can largely compensate the a em suppression, 
leading to an 0(1) correction to the decay amplitude to 
the 7 = 3/2 final state. This fact, however, does not spoil 
the possibility of extracting information on the UT with 
small hadronic uncertainty, as explained below. 

Dalitz plot analyses combined with isospin have al- 
ready shown their effectiveness in the extraction of a 



from B — > mm 9]. A proposal relying on isospin for ex- 
tracting 7 using a global analysis of B — > K s(7r7r)j= 0,2, 
including time-dependent CP asymmetries at fixed values 
of the Mandelstam variables, can be found in ref. 01 . We 
focus instead on K*n final states which allow us to per- 
form a simpler analysis with no need of time-dependent 
measurements. 



EXTRACTING CKM MATRIX ELEMENTS 
FROM B -» A'tttt DALITZ PLOTS 

Let us first illustrate our idea for the simplified case in 
which we neglect EWP contributions. To this aim, we 
write the amplitudes using isospin symmetry, in terms of 
Renormalization Group Invariant (RGI) complex param- 
eters jjjji obtaining 

A(K*+tt-) = V^VuPi-VZMEr-P?™) 
V2A(K*\°) = -V&Vuft-V&V^Ei + I?™ 1 ) 

V2A(K* + 7T°) = VtlVtsP! (1) 

-V: b V us (E 1 +E 2 + A 1 - Pf lM ) 
A(K*\+) = -V t * b V ts Pi + V: b V us (Ai-Pi lM ), 

where p^ GIM ^ represent (GIM-suppressed) penguin con- 
tributions, A\ the disconnected annihilation and E\ (E2) 
the connected (disconnected) emission topologies. B de- 
cay amplitudes are simply obtained by conjugating the 
CKM factors Vij . Similar expressions hold for higher K* 
resonances. 

Considering the two combinations of amplitudes 

A° = A(K* + TT-) + V2A(K* ir ) 

= -V: b V us (E 1+ E 2 ), (2) 

A = A(K*-n+) + V2A(K*°n a ) 

= -V ub V: s {E x +E 2 ), (3) 



the ratio 



R° 
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v: b v us 
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provides a clean determination of the weak phase 7. 
We now discuss how to extract Aq and Ao- Look- 
ing at the decay chains B° — > K* + (-^ K + tt°)tt" and 
B° —>■ K*°(-^ A'+7r~)7r°, one can obtain A from the 
_R" + 7T~7r° Dalitz plot, including the phase in a given con- 
vention, for example ImA(K* + ir~) = 0. Similarly, A° 
can be extracted from the K~ir + ir° Dalitz plot using 
the same procedure, choosing ImA(K*~ir + ) = 0. How- 
ever, in general, this choice does not reproduce the physi- 
cal phase difference between A(K* + tt~) and A(A"*~7r + ), 
which has to be fixed using additional information. 

This information can be provided by the Kg-K^ix^ 
Dalitz plot, considering the decay chain B° — > K* + (-^ 
K a ir+)Tr- and the CP conjugate B° -> 
K°ir~)ir + . These two decay channels do not interfere 
directly on the Dalitz plot, but they both interfere with 
the decays B,B — > p°(-^ ir + Tr~)Ks and with other res- 
onances contributing to the same Dalitz plot. There- 
fore the Dalitz analysis of B, B — > Kstt + tt~ should in- 
clude the p°Ks final state. In a time-integrated anal- 
ysis, the p a Ks final state comes from a mixture of 
B and B mesons, while the K* + ^-k~^ final state 
only originates from B (B) decay. Looking at the 
phases of A(K* +{ --^- (+ ^) relative to A(p°K s ), we can 
extract the phase difference between A(K* + ir~) and 
A(K*~ tt + ) with no need of resolving the flavor of the 
B in A(p°K s ) |3. 

A similar isospin relation involves charged B decays. 
We have 



(5) 
(6) 

(7) 



A+ = A(K*°ir+) + V2A(K* + ir°) 

= -V: b V ua {E 1+ E 2 ), 
A- = A(K*°tt-) + V2A(K*-n Q ) 

= -VubV;^ + E 2 ) , 

and the ratio 

A~ 
A+ 

As before, A^ can be extracted from the decay chains 
B ± K *±(_> K **)* and B ± -> K*°(-^ K ^ )^ 
entering the Kstt^tt Dalitz plot. Electric charge forbids 
the extraction of the relative phase of the two Dalitz plots 
along the way discussed above, so that a strategy based 
on theoretical arguments has to be adopted. In particu- 
lar, one can follow two possible paths. 

The first one is to use isospin symmetry to relate 
charged and neutral B decays: 

A{K*+tt-) + V2A(K*\°) 
-A(K*°TT + )-V2A(K* + Tr°) = (8) 
A(K*-n+) + V2A(K*°n°) 
-A(K*%-)-V2A(K*-7r ) = (9) 

so that the relative phases can be fixed. In this way, 
no additional information on 7 can be extracted from 



. However, the full information coming from charged 
and neutral -B's can be combined to improve the accu- 
racy of the determination of 7, thanks to the increase of 
available statistics. The second possibility is to use the 
penguin-dominated channel K*°tt + to fix the phase dif- 
ference between the amplitudes in the two Dalitz plots. 
In this way an independent, albeit more uncertain, de- 
termination of 7 can be obtained from RT . To illustrate 
this point, let us write down the phase for the A*°7r ± 
final state 



arg 



(A(K* tt+)) = /3 s +axg(l + ^|^A+e«A + 

V V tb y ts 



arg (A(K*\-)) = -ft + arg ( 1 + "g^A 



V ub V r 



US A + 2(5 



A + 



where ft = arg(-V ts V t l/ (V cs V; b )) and 

A x - Pf™ 



A+ e i5 A+ 



Pi 



(10) 



We now take advantage of the fact that 
|V^,K,s|/|V^Vy C 1 to simplify the above equa- 
tions and we obtain 



arg(A(A*°7r-)) = arg (A(K*°ir+)) - 2ft 

+2A+Im-^ 7 Jcos ( 5 A+ 
Vtb V ts 



(11) 



On general grounds, we expect A + ~ 0(1)- The error 
induced by the last term in eq. (fTTf) can be estimated at 
the level of \V ub V* s \ / \V tb V t * s \ ~ A 2 . The determination 
of 7 from R T is not as theoretically clean as the one 
obtained from R°. Nevertheless, the uncertainty induced 
by our dynamical assumption, being of 0(A 2 ), is much 
smaller than the expected experimental error (at least in 
the near future). 



INCLUSION OF ELECTROWEAK PENGUINS 

The inclusion of the effect of EWP's completely 
changes eqs. (0} and (7J . In fact, even though EPWs give 
a subdominant contribution to branching ratios (because 
of the 0{a em ) suppression with respect to the strong pen- 
guin contribution), they provide an 0(1) contribution to 
R° and R T (more generally, they provide an 0(1) cor- 
rection to CP violating effects in charmless b — * s de- 
cays). Fortunately, as we shall discuss in the following, 
the dominant EWP's (i.e. left-handed EWP operators) 
can be eliminated at the operator level, so that no addi- 
tional hadronic matrix elements are introduced. The net 
effect of EWP's is that R° and R^ depend not only on 
7 but also on other CKM parameters. 

Let us consider the effective Hamiltonian for b — ► s 
transitions given for instance in eq. (5) of ref. There 
is a hierarchy in the values of the Wilson coefficients for 
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FIG. 1: Bounds in the p-rj plane from the analysis of B° decays, assuming a measurement of the relative phase with an error 
of 20° (left) or 40° (right). The output of the present UT fit analysis [l| is shown as a reference. 
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EWP operators: |Cg,io| ^ |C7,s I - Let us therefore ne- 
glect the effect of Qi^ and focus on Qg^o, barring a very 
large dynamical enhancement of {Qr$) in the case of B 
decays. There is an exact operator relation that allows 
to eliminate Qg t io Il2l in favour of current-current and 
penguin operators [17| : 



Q9 — 2 (Q™ 



ho = g(Q 



suu 

1 



Q s 2 cc )+m cc - 
q? c ) + ?>q? c -\qi 



(12) 



so that the effective Hamiltonian becomes 



71 



v: h v us c+ - -v; b v ts c™ ) (q s a 



-Q 



scc\ 

+ I 



v: b v tts c_ + -v; b v ts c™ ) - 



v; b v ts H AI =° 



(13) 



where Q± = (Qi ± Q 2 )/2, C± 
Cg ± C10 • We observe that the relation 



Ci ± C 2 and CJ W 



C2 



c 4 



C- 



(14) 



is exact at the LO and is broken by 0{a s a em log) cor- 
rections. Numerically, eq. (|14fl holds with high accuracy, 
being violated at the percent level. Using this relation, 
we can write 



V: b V us (l + K EW ) C+(Q 



Q 



scc\ 

+ ) 



-C- (Ql" u - Q s 



- VXV ts H 



A/=0 



,(15) 



where 



3 C ew 
~2~C^ 

3C ,EW 

2~c7 



v t * b v ts 



ub 



1-A 2 



A 2 (p + irj) 



0{\ 2 



Therefore R° and R T can be rewritten as 



1 Ktt.\ 



R ^ 
where 



-2i(7+arg(l+K EW )) _ 



1 + K* 



EW ° + 



1-kew C- rc i6~ 
1+kew C+ 



(K*ir(I = 3/2)\Q_\B) 
(K*ir(I = 3/2)\Q+\B) 



(16) 



(17) 



(18) 



While for B — > Ktt decays the SU (3) flavour symm etry 
guarantees that (Ktt(I = 3/2)|Q_|i?) vanishes pi Il3l|. 
the same argument, based on the symmetry property of 
the final state wave function, does not apply to K*n final 
states |l4j . However, using factorized amplitudes and 
form factors as given in ref. 0, one obtains 



f Z?-B-»7T f aB—>K' 



f K >F B ^ + UA*- 



■K~ 



< 0.05. 



(19) 



While this numerical result depends on the estimate of 
the form factors, the good agreement between QCD sum 
rules and lattice QCD calculations makes it rather ro- 
bust. In our analysis, however, we do not assume a spe- 
cific model for computing the amplitudes, rather we let 
r to vary in the conservative range 0-0.3. 
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Concerning k ew , using the values C + (™b) = 0.877, 
Cl w {m b ) = -1.017 a em and p = 0.216, fj = 0.342 Q, we 
obtain k ew = —0.35 + 0.53 i. One can thus verify that 
k ew is an 0(1) correction to the decay amplitude to the 
7 = 3/2 final state. 

The above equations allow to translate the experimen- 
tal results for R° and R T into allowed regions in the p-fj 
plane. Neglecting terms of 0(r 2 ), for a given value of R° 
one obtains the linear relation fj = — tan(-| argi?°)(p — 
po), where 

3Cf w 12Cf w C^r cos r ~ 

P ° ~ ~ |_2C++3CT " (2C+ + 3Cr) 2 . ' 

i^ + G(A 2 ) (20) 

Including terms of 0(r 2 ) or higher, the relation between 
fj and p becomes quadratic, but the deviation from the 
result in eq. (121 III is small and mainly amounts to a 
depletion of the region fj ~ 0. Furthermore, the ef- 
fect of r is constrained by the measurement of |i? 0,T |, 
since|i?° :=F | — 1 oc rsin# r . 

We can test this new idea using the experimental re- 
sult of the BaBar Collaboration on B° — > K + ir~ir (> de- 
cays [l^. Among the other measurements, this analysis 
provides the decay amplitudes for B° and B° decays to 
K*n and to A"*(1430)7r final states. We implement this 
information in our analysis using directly the shape of 
the multidimensional likelihood from BaBar (including 
all correlations). In this way, we obtain an error of 38° 
on argi? . Unfortunately, the Kstt + tt~ Dalitz plot is not 
yet available so that we cannot fix the relative phase of 
B and B decays at present. For the sake of illustration, 
we assume a central value for this relative phase such 
that the constraint on p and fj from eq. (fTTjl is compati- 
ble with the SM UT fit result Q. For the experimental 
uncertainty on the relative phase, we consider two cases, 
corresponding to ±20° or ±40°, leading to the constraints 
exhibited in Fig. ^ obtained without expanding in r. 

The situation can be improved by fitting R.%*^ and 
-fi^*(i430)7r directly from data, cancelling out part of the 
systematic error. In addition, the determination of the 
UT parameters can be further improved with the experi- 
mental measurement of Rk,^ and ^J*( 1 4 3 o) Tr an d adding 
Belle data. 

Let us finally comment on the sensitivity to NP of our 
analysis. Making the very reasonable assumption that 
NP effects only enter at the loop level, we can envisage 
three possibilities. First of all, NP could affect the co- 
efficients of QCD penguin operators. In this case, the 
analysis of R° is completely unaffected, while the phase 
of NP contributions would modify eq. 1|11|) . This could 
produce a discrepancy between the constraints on the UT 
obtained from R° and R T using eq. ifTTJ) . A second pos- 
sibility is that NP modifies EWP coefficients, respecting 
however the hierarchy Cq^q 3> CV,s- I n this case, the only 



effect would be a modification of k ew , so that the con- 
straint on the UT obtained using the SM value for k ew 
could be inconsistent with the SM UT fit result. Finally, 
NP could produce contributions to EWP operators such 
that Cg^o ~ CV,8, or give rise to new Al — 3/2 opera- 
tors that cannot be eliminated. In this case, one would 
observe \R°^\ ^ 1. Present data give = 0.96 ±0.17. 
A small |i?°' T | — 1 could also be generated by (Q-) or by 
a large dynamical enhancement of (Qt,s) within the SM. 

CONCLUSIONS AND OUTLOOK 

We have presented a new method to constrain the Uni- 
tarity Triangle using B — > Kmr decays. This can be 
achieved with both neutral and charged B decays, using 
the amplitude ratios R° and R T defined in eqs. Q and 
(JZJ) . The theoretical uncertainty is negligible with respect 
to the foreseen experimental error. We have discussed 
in detail how to take into account electroweak penguins. 
Our exploratory study shows that this new constraint can 
have a sizable impact on the Unitarity Triangle analysis 
in the near future. We have discussed possible improve- 
ments and sensitivity to New Physics. 
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